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This  paper  considers  the  problem  of  maximizing  the  energy  or 
average  power  transfer  from  a  nonlinear  dynamic  n-port  source. 
The  main  theorem  includes  as  special  cases  the  standard 
linear  result  Yioacj  =  Y*source  and  a  recent  finding  for  non¬ 
linear  resistive  networks.  An  operator  equation  for  the  optimal 
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I.  Introduction 


This  paper  addresses  the  problem  of  extracting  the  maximum  energy  or 

1 

average  power  from  a  source  with  the  topology  shown  in  Fig.  1.  As  in 
[1],  the  problem  is  formulated  as  finding  the  optimal  output  voltage 
v(-)  for  each  current  source  waveform  i  (•)  rather  than  finding  a  load 
that  maximizes  the  power. 

The  central  result  is  the  operator  equation  (6)  for  v(*)>  Theorem  1 
gives  conditions  that  guarantee  uniqueness  and  global  optimality  of  the 
solution:  the  standard  result  for  linear  systems  [1]  and  recent  work 
on  resistive  nonlinear  systems  [2]  follow  as  special  cases.  Equation  (11) 
defines  a  practical  algorithm  for  solving  (6),  and  Theorem  2  gives  conditions 
that  guarantee  convergence. 

The  solution  v(»)  can  be  of  engineering  value  in  two  ways.  First,  the 
average  power  P(v)  tells  us  the  optimal  performance  that  is  possible  in 
principle.  Second,  v(*)  itself  is  a  concrete  design  goal.  If  the  source 
admittance  operator  F  is  continuous,  a  load  for  which  the  output 
approximates  v(«)  (in  the  Hilbert  space  norm  used  in  this  work)  will  absorb 
an  average  power  that  approximates  F(v). 


1.  Reference. [1]  actually  deals  with  the  dual  network,  where  the  source 
appears  in  Thevenin  form. 
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II.  Results 


2.1)  Notation  and  Definitions 

Let  L  be  any  real  inner  product  space  and  L  any  linear  subspace  of  L. 

An  operator  F:  L  +  L  is  said  to  be 

a )  strictly  increasing  i f 

<F(y)  -  F(x) ,  y-x>  >  0,  Vx  ?  y  e  Ls  (1) 

b)  uniformly  increasing  if  for  some  6  >  0, 

^F(y)  -  F(x),  y-x)  >  6||y-x||2,  Vx,y  e  L,  (2) 

c)  Lipschitz  continuous  if  for  some  K  •>  0, 

II  F(y)  -  F(x)  ||  <_  K||  y-x||  ,  Vx,y  e  L.  (3) 

Let  Ls  L1  be  any  real  inner  product  spaces  and  L(L,L')  denote  the  space 
of  continuous  linear  maps  from  L  to  L',  with  the  operator  norm  [3,  p. 53] . 
For  AeL(L ,L 1 ) s  let  Aac^  denote  the  adjoint  of  A. 

Given  an  operator  F:  L  L1  and  x,heL,  suppose  there  exists  an  element 
denoted  6F(x,h)  of  L1  such  that 


1  im 
t+0+ 


F(x+th)  -  F(x)  -  6F(x,h) 
t 


L' 


=  0. 


Then  6F(x,h)  is  called  the  Gateaux  variation  of  F  at  x  for  the  increment  h 
[4 , p . 251  ] .  If  SF(x,h)  exists  for  all  x,heL,and  if  for  each  xeLthe  map 
h  ^  5F(x,h)  is  an  element  of  L(L,L'),  then  F  is  said  to  be  Gateaux 
differentiable  on  L.  In  this  case  the  map  x  ->  6F(x,*)  is  called  the  Gateaux 


derivative  of  F  and  denoted  DF:  L  -*  L(L,L')  [4, pp. 255-256],  Similarly 
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6 F ( x , • )  is  denoted  DF(x)eL(L,L‘ ) ,  and  <5F(x,h)  is  denoted  (DF(x) ) h e  L 1 .  The 

value  of  using  the  Gateaux  derivative  rather  than  the  more  restrictive 

Frechet  derivative  [4, Chap. 3]  will  become  apparent  in  section  3.1. 

The  Hilbert  space  L^  is  the  set  of  all  measurable  functions  x:  IR-*  JRn 

o 

such  that  the  integral  of  x.(*)  over  IR  is  finite,  j=l,,..,n,  equipped 

J 

A  1  /9 

with  the  usual  inner  product  and  norm,  ||x||  ^=(x,x) '  . 

2 

For  each  T  >0,  y  is  the  set  of  all  periodic  measurable  functions 
x:  IR  -*  3R  n  with  period  T  such  that  the  integral  over  one  period  of  x.(*)  is 
finite,  j=l,...,n.  It  is  a  Hilbert  space  with  the  "average  power"  inner 
product 


fT 

x(t) *y(t)  dt. 


0 


where  x*y  is  the  Euclidean  inner  product  on  ]Rn. 
II  x||  T  4(x,x)](2. 


The  norm  on  L 


2 

n  ,T 


(4) 

is  denoted 


2.2)  Main  Theorem 

Theorem  1  (Maximum  Average  Power  in  the  Periodic  Steady  State) 

Fix  T  >  0  and  let  the  n-port  M,  in  Fig.  1  be  characterized  by  an 

2  2 

admittance  operator  F:  Ly  -*•  Ly,  where  Ly  is  any  linear  subspace  of  Ln  y. 
Suppose  F  is  Gateaux  differentiable  on  Ly  and  that  the  associated  operator 
H:  Ly  ->  Ly,  given  by 

adj 

H:  v  h-  F(v)  +  (DF^vj)  v  ,  (5) 

.  ^  *  2 

2.  Thus,  if  y(*)  has  period  T  and  lies  in  L  T,  the  response  i(*)  of  W, 
cannot  have  subharmonics.  ’ 
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is  strictly  increasing. 

/s 

Then  for  each  H(Ly)  there  is  a  unique  solution  v(is)eLy  to 

is=H(v),  (6) 

3 

and  the  average  power  absorbed  by  the  load* 

P(v)  -{i0,v)T  =  (i$  -  F(v),v)y  »  (7) 

has  a  unique  global  maximum  over  Ly,  which  is  attained  at  v  =  v(is). 

Gorol lary  (Maximum  Total  Energy  for  Transients) 

2 

Let  L  be  a  linear  subspace  of  Ln,  and  substitute  L  for  Ly  in  the  assumptions 

&  - 

of  Theorem  1.  Then  the  same  conclusions'  hold,  but  with  v(i  )eL  maximizing 

the  total  energy  E(v)  ^  -  F(v),v^  over  L. 

Note  that  in  general  F  can  be  nonlinear  and  time- varying. 

In  applications  one  might  wish  to  restrict  attention  to  currents  and 
2 

voltages  in  Ln  T  with  additional  properties  such  as  continuity  or  boundedness, 
n  j  I 

o 

This  is  the  reason  for  introducing  LjC  Ln  y  in  the  formulation  of  Theorem  1. 

The  essential  idea  behind  the  theorem  is  that  a  solution  v(*)  of  (6) 
is  a  stationary  point  of  F:  Ly  -»■  IR  ,  and  the  monotonicity  assumption  on  H 
guarantees  that  F  is  strictly  concave.  Details  follow. 


3.  A  more  explicit,  but  cumbersome,  notation  would  be  P ( v , | s ) -  Using  it. 
Theorem  1  states  that  Vy,iseLj,  P(v,is)  <  P(v(is),is)  if  v  t  v(is). 

4.  For  the  Corollary,  the  adjoint  is  of  course  taken  with  respect  to  the 
inner  product  on  L^  rather  than ,  ^y. 
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Proof  of  Theorem  1 

Uniqueness  of  the  solution  to  (6)  follows  from  the  fact  that  H  is 
strictly  increasing.  By  the  chain  rule  for  the  composition  of  Frechet- 
and  Gateaux-differentiable  functions  [4, p.253]  F  is  Gateaux  differentiable 
and  for  all  x,he  Ly, 

(DP(x))h  =^is  -  F(x),h^y  -  ((DF(x))h,x^y  = 

(is  ■  -  (D!(?))adj  ?’H)t  = 

<1S  -  H ( x ) , h)T .  (8) 

Thus  if  is£H(LT), 

a)  Dr(!(is)}  =?£L(lt>lt)> 

b)  given  any  x,y  eLy,the  map  x  -*■  F[x+x(y-x)]  is  differentiable  at 
each  XeIR  s  and 

c)  ^  P[x+A(y-x)]  =  ^is  -  H[x+x(y-x)] ,  (y-x))y. 

To  show  that  v ( i s )  globally  optimizes  F,  fix  is£H(Ly),  let  v  =  v(i‘s), 
and  choose  anyVELy,  v  f  v.  Then 

P(v)  -  P(v)  = 

F[v+x(v-v)]j  -  F[v+x(v-v)] 

~  ~  I X=1  ~  ~  ~  x=0 

rl 

P[v+X(v-v)]}  dx.  (9) 

Using  c),  the  integrand  above  is 
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^is  -  H[v+a(v-v) ] ,  v-v)T  = 

(since  i$  =  H(v)) 

^^Htv+xlv-v)]  -  H(v) ,  [v+x(v-v)]  -  [v]^T,  Va> 0, 

and  the  integrand  vanishes  at  A=0.  The  inner  product  above  is  strictly 
positive  for  \?0  since  H  is  strictly  increasing  by  assumption.  Thus  the 
integrand  in  (9)  is  negative  for  A>0  and  zero  for  a=0,  so  P(v)  <  F(v)  as 
claimed.  £ 

The  proof  of  the  Corollary  is  essentially  identical  and  will  be  omitted. 

2.3)  Relation  to  "Impedance  Matching"  Ideas 

The  emphasis  in  this  paper  is  on  finding  the  optimal  output  voltage 
v(*)»  not  the  optimal  load.  But  the  relation  to  impedance  matching  ideas 
deserves  comment. 

If  the  load  in  Fig.  1  is  taken  to  be  the  (generally  noncausal)  admittance 
Gopt:  LT  LT’  defined  by 

adj 

Gopt:  v  *  (DF(v))  v  ,  (10) 

then  the  network  is  uniquely  solvable  given  any  i  e  H(Ly),  and  the  output  voltage 
v(*)»  which  necessarily  equals  v(is),  globally  optimizes  F.  This  generally 
noncausal  load  is  "matched"  to  the  source  for  al  1  inputs  i  eH(Ly),  and  this 
result  holds  generally  for  a  nonlinear,  time-varying,  even  noncausal  source 
admittance  F.  The  reader  can  easily  verify  that  in  the  LTI  case  (10)  reduces 
to  the  standard  linear  theorem  Y]oac|(j“)  =  ^source^^’  ^ore  deta^  f°r 
the  linear  1-port  case  is  given  in  Section  3.1. 
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Of  course  in  practice  one  has  a  causal  load,  usually  predetermined,  and 
wishes  to  couple  it  to  the  source  through  a  lossless  matching  network  designed 
to  maximize  the  absorbed  power  over  a  range  of  inputs.  In  the  linear  case 
this  important  problem  is  called  "broadband  matching"  [5-8].  We  note  that 
in  both  the  linear  and  nonlinear  cases  the  problem  can  be  viewed  as 
compensating  or  coupling  to  a  predetermined  load  using  lossless  elements  in 
such  a  way  that  the  response  approximates  that  of  the  noncausal  exact  match 
Gopt  over  tfie  1nput  ran9e  interest. 

For  a  particular  drive  is,  the  situation  is  somewhat  different.  The 
optimal  voltage  v(*)  is  unique,  but  the  optimal  load  is  not:  the  only 
requirement  on  G  is  that  G(v)  =  Gopt(v).  In  the  linear  case  where  F  and  G 
are  respectively  represented  by  admittance  matrices  Y  (jw)  and  Y^(ja)),  there 
are  in  general  infinitely  many  optimal,  positive  semidefinite  choices  of  YL 
at  a  given  w  for  which  the  network  is  uniquely  solvable  [9].  The  problem  of 
finding  solutions  in  particular  classes,  such  as  the  class  of  resistive 
loads,  is  studied  in  [10]. 

2 .4  )  Numerical  Algorithm 

Equation  (8)  shows  that  i  -  H(v)  is  the  gradient  [3,p.72],  [4,p.l96] 
of  P  at  v,  V  i$,v  e  Lj.  This  suggests  that  we  attempt  to  maximize  F  by  a 
simple  "hill-climbing"  algorithm  of  the  form 

Xj+1  =  A(is  -  H ( x j ) )  +  x.  =M(Xj)  (11) 

for  some  a>0.  Note  that  under  the  assumptions  of  Theorem  1,  vf  x.  -*•  x  £  Ly 
and  H  is  continuous,  then  i  =  H(x)  and  x  globally  maximizes  P.  By  tightening 
the  assumptions  a  little  further,  we  can  guarantee  convergence  for  all 
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sufficiently  small  positive  A . 

Theorem  2 

Strengthen  the  assumptions  of  Theorem  1  by  supposing  further  that  Ly 

is  closed  and  H  is  uniformly  increasing  and  Lipschitz  continuous  on  Ly. 

(See  (2),  (3).)  Then  for  any  i$£Ly,  any  initial  guess  xQeLy,  and  any 
2 

Ae  (0,  26/K  ),  the  sequence  generated  by  (11)  converges  to  v(is). 

Remark 

Note  that  Theorem  2  also  guarantees  existence  of  a  solution  to  (6)  for 

all  iseLy,i.e. ,  H(Ly)  =  Ly. 

Proof 

2  2 

Since  LyC  Ln  y  is  closed  and  l_n  T  is  complete,  Ly  is  complete  [11]. 

It  remains  to  show  that  M  is  contractive,  i.e.,  that  for  some  C  <  1, 

II  M(y)  -  M(x)  1 1  y  <  C  ||  y-x  1 1  y,  Vx,y  e  Ly,  (12) 

to  guarantee  II  xn  -  v(i  ) 1 1  y  ->  0  by  the  contraction  mapping  theorem  [3,p.l02], 
[12,p.28].  But 

II  M(y)  -  M(x)  ||  2  = 

~  ~  y 

(y  -  x  -  A(H(y)  -  H(x) ) ,  y-x  -  A(H(y)  -  H(x)))y  = 

||  y  -  x|  |  2  -  2A(H(y)  -  H(x),  y  -  x^)y  +A2||  H(y)  -  H(x)  ||  2  < 
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and 


(1  -  2X6  +  A2)  ||  y-x ||  2  =C2(x)  1 1  y-x|  |  2  , 

~  ~  T  ~  ~  T 

C2(x)  <  1,  VXe(0,26/K2). 


III.  Examples 


3.1)  Linear  Operators  and  Memoryless  Operators 

2 

Consider  the  time-invariant  scalar  case  for  simplicity,  and  let  LT  stand 
2 

for  L,>r 

If  F  is  the  convolution  operator:  vr  a*b  where  a:  F  +  1R  is  absolutely 
a  2  2 
integrable,  then  for  each  T  >  0,  F£  is  a  continuous  linear  operator:  Ly  -»■  Ly 

A 

and  therefore  Gateaux  (in  fact,  Frechet)  differentiable.  Since  F  is  linear 

X/ 

DFn(x)  =  F  ,  and  the  reader  can  easily  verify  that  (DF  (x))adj  = 

-s « 4  ** 

aaj 

F  :  v(-)  h-  a(-*)*v-(*)»  i.e.,the  adjoint  operation  turns  the  impulse  response 

X/ 

around  in  time.  Furthermore,  H  :  v(*)  h-  [a(  * )  +  a(- » ]Q*  v  ( • )  is  strictly 
increasing  on  Ly  for  each  T  >  0  iff  Re  {a(ju>)}  >  0  for  all  w,  where  a  is  the 
Fourier  transform  of  a.  This  follows  from  a  slight  modification  of  [12:pp.25, 
174,235].  Similar  results  hold  if  a(*)  contains  impulse  functions  as  well 
[12:pp. 246-247].  Thus  v(*)  h-  a(-*)*v(*)»  and  Gopt  is  represented  in  the 

frequency  domain  by  the  complex  admittance  a  (jw).  Therefore  Theorem  1  and 
equation  (10)  reduce  to  the  standard  result  Yi oacj ( J “)  =  Ysource^°^ 
is  linear  and  time-invariant. 

Suppose  Fm  is  memoryless  but  possibly  nonlinear,  i.e.,  N-j  is  a  resistor 

with  the  constitutive  relation  i=b(v).  Assume  that  b:  K -*  IR  is  differentiable 

and  its  derivative  b'(*)  is  bounded.  Then  b  is  Lipschitz  continuous  on  F 

2  2 

and  hence  for  each  T  >  0  the  operator  Fm:  v(t)  h-  b(v(t))  maps  Ly  into  Ly. 

Using  Prop.  13  of  [13: p. 85]  and  the  Lebesgue  Convergence  Theorem  [13:p.88], 
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"  2  2 
one  can  show  that  Fm  is  Gateaux  differentiable  on  Ly  and  that  for  all  x,y  e  Ly, 

(DFm(x))y  =  y(*)b' (x(*))eLy-  Furthermore  if  h:  v  h-  b(v)+vb'(v)  is  a 

strictly  increasing  function  on  1R  ,  then  Hm:  v(t)  h-  b(v(t)')  +  v(t)b'(v(t)) 

2  2 

is  a  strictly  increasing  operator:  Ly  ->  Ly.  Thus  Theorem  1  reduces  in  this 
case  to  the  result  in  [2]. 

2  2  2 

The  reader  can  easily  check  that  DFm:  L^:  +  L(L^,L^)  is  not  continuous 

unless  b' ( *)  is  constant.  Thus  if  W-,  is  a  resistor  with  any  nonlinearity 

(other  than  the  trivial  i=gv  +i),  Fm  is  not  Frechet  differentiable  [4, Chap. 3] 

2 

on  Ly.  This  is  the  reason  Theorem  1  was  formulated  in  terms  of  the  weaker 
Gateaux  derivative. 

3.2)  Positive  Linear  Combinations  of  Operators 

The  (noncausal)  matched  load  (10)  for  the  source  admittance  F  is  related 

adj 

to  F  by  a  mapping  £,  £(F)  =  Gopt:  vh-  (DF(v))  v  .  Note  that  £  is  linear; 
i.e.  A(aFy  +  bF2)  =  a  £(F-j)  +  b£(F2).  Given  F-j  and  F2:  Ly  ->  Ly,  consider 
F  =  aF-|  +  bF2>  The  reader  can  easily  verify  that  if  F^ ,  F2  satisfy  the 
conditions  of  Theorem  1  (resp.  Theorem  2),  then  F  also  satisfies  Theorem  1 
(resp.  Theorem  2),  provided  a  >_  0,  b  0,  a  +  b  >  0. 

For  example,  consider  the  source  shown  in  Fig.  2,  where  W-j  consists  of  the 
parallel  connection  of  an  LTI  1-port  and  a  nonlinear  resistor.  If  Y  and  g 
satisfy  the  conditions  in  section  3.1),  then  the  (noncausal)  matched  load 
has  the  form  shown  in  Fig.  2. 

3.3)  Circuit  Example 

Suppose  the  source  takes  the  specific  form  in  Fig.  3,  with  the  resistor 
curves  shown  in  Fig.  4.  The  convolution  kernel  a(t)  =  e-t,  t  >  0,  for  the 
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series  connection  of  inductor  and  resistor  satisfies  the  assumptions  of 
section  3.1.  The  resistor  curves  g^  are  differentiable  everywhere  and 

hk(v)  ^gk(v)  +  vgk'(v)  =  ( k+1 ) v | v | k_1 ,  k  =  1,2,3,  (13) 

with  h-j(O)  =  0.  All  the  assumptions  of  section  3.1  are  satisfied  except 
that  the  derivatives  g2 1 ( * )  and  g3 1 ( • )  are  unbounded.  (Since  they  are  bounded 
on  every  bounded  subset  of  1R  ,  a  more  detailed  argument,  omitted  here, 

—  oo  ? 

shows  that  the  solutions  obtained  below  maximize  P  over  l_T  0  Ly,  which  is 
certainly  sufficient  in  practice.) 

To  find  the  optimal  output  v  in  the  three  cases,  we  carried  out  the 
iterative  procedure  (11),  which  becomes  in  this  instance 

Vi(t)  =  ■ 


6  sin(t)  -  (k+l)x.(t) |x.(t) 

J  J 


k-1 


e  ^”Tb.(T)dT 

J 


+  Xj(t),  k  =  1 ,2,3. 


L  •'-«  J  (14) 

Miss  Pearl  Yew  of  MIT  has  written  a  program  in  PASCAL  to  do  the  numerical 
solution.  It  was  run  on  the  DEC 20  in  MIT’s  Research  Laboratory  of  Electronics 
with  an  initial  guess  of  xQ(*)  =  0,  and  found  to  converge  fairly  rapidly 
for  small  positive  values  of  X.  The  results  are  shown  in  Fig.  5. 

Since  g-|  represents  a  linear  resistor,  it  follows  from  the  traditional 
linear  theorem  that  v(t)  =  2sin(t)  for  k=l ,  in  agreement  with  the  numerical 
solution.  Note  that  the  instantaneous  current  drained  by  the  nonlinear  source 
resistor  increases  in  magnitude  with  k  for  |v|  >  1  but  decreases  for 
| v |  <  1.  Thus  it  is  intuitively  reasonable  that  the  optimal  output  spends  a 
progressively  greater  percentage  of  time  in  the  region  |v|  <  1  as  k  increases. 
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Figure  Captions 


Fig.  1  The  solution  of  the  operator  eguation  (6),  given  a  particular  i  (•)» 
is  the  optimal  output  voltage  v(*)-  It  can  be  achieved 
with  a  variety  of  loads. 

Fig.  2  The  optimal  load  admittance  is  obtained  by  a  linear  operator  £ 
on  the  source  admittance.  Thus  the  optimal  load  for  a  paralfel 
connection  of  source  admittances  is  the  parallel  connection  of  the 
optimal  loads  for  each  source  separately. 

Fig.  3  Theorems^!  and  2  let  us  numerically  determine  the  optimal  output 

voltage  v(*)  for  this  circuit  when  the  resistor  curves  are  as  shown 
in  Fig.  4. 

Fig.  4  The  three  resistor  curves  for  the  circuit  in  Fig.  3  are  g,  (v)  -  vlvl^, 
k=l  ,2,3,  with  3,(0)  iO.  k 

Fig.  5  One  period  of  the  optimal  output  voltages  for  the  circuit  in  Fig.  3. 
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